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In this paper, we continue to study a unified dark fluid model with a constant adiabatic sound 
speed but with the entropic perturbations. When the entropic perturbations are included, an 
effective sound speed, which reduces to the adiabatic sound speed when the entropic perturbations 
are zero, has to be specified as an additional free model parameter. Due to the relations between 
the adiabatic sound speed and equations of state (EoS) (?s^ad{o) = wia) — dln(l + ui(a))/3dlna, the 
equation of state can be determined up to an integration constant in principle when an adiabatic 
sound speed is given. Then there are two degrees of freedom to describe the linear perturbations 
for a fluid. Its micro-scale properties are characterized by its EoS or adiabatic sound speed and an 
effective sound speed. We take the effective sound speed and adiabatic sound speed as free model 
parameters and then use the currently available cosmic observational data sets, which include type 
la supernova Union 2.1, baryon acoustic oscillation and WMAP 7-year data of cosmic background 
radiation, to constrain the possible entropic perturbations and the adiabatic sound speed via the 
Markov Chain Monte Carlo method. The results show that the cosmic observations favor a small 
effective sound speed cl^^ff = 0.00155tE!;^^5!.-|)°oS?5f--"o°0M5f in l,2,3cr regions. It means that a 
UDF model with small entropy perturbation is favored but the pure adiabatic case is not ruled out. 



I. INTRODUCTION 

In the last few years, the unified dark fluid models 
p]-[l^ were investigated as a possible explanation to an 
accelerated expansion phase of our Universe [H, [l4|. 
These models are inspired by the facts that above 96% 
of the energy content in the Universe is made of un- 
known dark component. These unified dark fluid mod- 
els include the popular generalized Chaplygin gas (gCg) 
model as a sample which is a generalization of the 
Chaplygin gas (Cg) model or a coined model from the 
ACDM model 0]. Actually, the EoS' of these unified 
dark fluid are specified in different models. Then their 
linear perturbation evolutions are determined via the re- 
lations ^^(a) = w{a) — dln(l -I- w(a)) /idhia between 
the adiabatic sound speed and its EoS of UDF when the 
entropic perturbations are zero. On the contrary, when 
the adiabatic sound speed is preposed, the EoS can also 
be determined up to an integration constant in princi- 
ple. Based on this point, the so-called AaCDM or CASS 
model was proposed 0, 0, US] where a constant adiabatic 
sound speed = a was assumed . The case of a = 
was discussed in Ref. [ll|. And the time variable sound 
speed cases were also discussed in [l^] . As pointed out by 
the authors [l^, when a is negative, this model can be 
seen as the attractor solution of quintessence scalar field 
dynamics. Also, it can be treated as k-essence scalar field 
Q. Interestingly, it can avoid the so-called averaging 
problem [iTj when the perturbations become nonlinear 
[3|. So, in this paper, we are going to investigate this 
unified dark fiuid model. 
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However when the entropic perturbations are included 
the adiabatic sound speed is not enough to characterize 
the micro-scale properties, one has to introduce the ef- 
fective sound speed as a specified model parameter as an 
addition to its EoS. For a generalized dark matter, the 
effective sound speed was defined in Ref. [ISj . The effec- 
tive sound speed is reduced to the adiabatic sound speed 
when the entropic perturbations vanish. 

The model with a constant adiabatic sound speed and 
entropic perturbations was investigated in Ref. ^is!], 
where the limited cases of effective sound speed c^j^ = 
a, 0, 1 were discussed. However, the effective sound speed 
should be a free model parameter to be determined by 
the cosmic observations instead of being fixed to a spe- 
cial value by hand. So, in this paper, we will consider a 
more general case where the effective sound speed 
is taken as a free model parameter in the range of [0, 1]. 
And we will use the currently available cosmic observa- 
tional data sets, which include type la supernova Union 
2.1, baryon acoustic oscillation and WMAP 7-year CMB, 
to determine the model parameter space via the Markov 
Chain Monte Carlo (MCMC) method. 

This paper is structured as follows. At first, in sec- 
tion |TT1 we give a brief review of the unified dark fluid 
model with a constant adiabatic sound speed (CASS) 

jjjj = a. In this section, the energy density and back- 
ground evolution equation are shown. In section iHll by 
using Markov Chain Monte Carlo (MCMC) method with 
currently available cosmic observational data sets which 
include type la supernova Union 2.1, baryon acoustic os- 
cillation and WMAP 7-year CMB, we show the parame- 
ter space. A summary is presented in section ITVl 
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II. A BRIEF REVIEW OF CASS MODEL 



following forms 



In this section, we will give a brief review of CASS 
model which has constant adiabatic sound speed — a, 
for the details please see the Ref. [1] . The energy density 
and equation of state (EoS) of this UDF are given in the 

I 



Pd = PdO 



(l-i?,) + i?,a-3(i+") 
(l-S,) + S,a-3(i+")' 



(1) 
(2) 



where Bs in the range < Bs < 1 and a are model 
parameters. In the pure adiabatic perturbation case, the 
value of a should be fixed in the range of [0, 1]. When the 
entropic perturbations are included, the adiabatic sound 
speed can be negative [isj . So, in this paper, we assume 
it is in the range [—1, 1]. When Bs is zero, this UDF is 
a cosmological constant. If a = is respected, this UDF 
is standard ACDM model at the background level. And 
when Bs — 1 and a — are set, it is cold dark matter. In 
a spatially flat FRW Universe, the Friedmann equation 
is given as 



= (1 - 1^6 - (1 - Bs) + Bsa 



'3(1+q) 



-3 



(3) 



where H is the Hubble parameter with its current value 
Hq = 100/ikm s~^Mpc~^, and fli (i = 5, r) are dimen- 
sionless energy parameters of baryon, radiation density 
respectively. 

In Ref. [1], the CASS model as a UDF (with- 
out any decomposition) was tested by using SN la, 
BAO and full information from WMAP-7yr data sets, 
where the model parameter space was obtained. Here 
we just show the 1,2a regions for a and Bs'. a — 

n nnn/1S7+0 000117+0. 000728 r> _ n OOQ+0-0134+0.0274 

U.UUU^O ( _0. 000487-0. 000487 '^'^'^^ — U.ZZa_o 0133-0.0249- 

Considering the perturbation in the synchronous 
gauge, the perturbed metric reads 



a^ir) [-dr^ + {5^j + /ly (x, T)dx'-dx^)\ , (4) 



where r is the conformal time and hij is the metric per- 
turbation. From the conservation of energy-momentum 
tensor Tj^.^ = 0, one has the perturbation equations of 
density contrast and velocity divergence for dark fluid in 
the synchronous gauge 

5d = -(1 + wd){ed + ^)- - wd)5d, (5) 

2 opd 
^ - 3cL.) + ^-^^kHd k^<yd (6) 

following the notations of Ma and Bertschinger [l9| . 
where the definition of the adiabatic sound speed 



•'s^ad 



Pd 
Pd 



Wd 



Wd 



'm{i + wd) 



(7) 
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is used. For the gauge ready formalism about the pertur- 
bation theory, please see [20]. For a pure barotropic fluid, 
it has an imaginary adiabatic sound speed which causes 
instability of the perturbations when its EoS is negative, 
for example the w = constant quintessence dark energy 
model. The way to overcome this problem is to allow 
an entropy perturbation and to assume a positive or null 
effective speed of sound. Following the formalism for a 
generalized dark matter 15], one can separate out the 
non adiabatic stress or entropy perturbation for the UDF 



PdTd = 5pd 



l,ad^Pd 



(8) 



which is gauge independent. In the rest frame of UDF 
by introducing the effective speed of sound eff^ the 
entropy perturbation is specified as 



wdrd = icl,fj-cl,^)sr'- 



The gauge transformation into an arbitrary gauge 



(9) 



- Sd + 3H{1 + Wd) 



fc2 



(10) 



gives a gauge-invariant form for the entropy perturba- 
tions. By using the Eqs ([8j) ([9]) and (|TOl) . one can recast 
Eqs. ©, and © into 
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6d = -(1 + Wd)i9d +1^)- m^lad - y^d)6d - mcleff - cladWd 



-(1 + wd){9d + -) + — ^<5, ~ m{cl,jf - cLJ'^r 



Wd 

-(1 + WdXOd + b + T^Sd - mcleff - cl^d) 

Z 1 + Wd 

^2 



5d + m{i + wd)-^ 



-nOd 



c 



l+Wd 



k^Ud 



l+Wd 



Sd + 3n{l + Wd)^ 



k^dd 



= -n{l - 3c; 



^eSd - k^ad 



1 + Wd 



please see also in [l^ and [2l| for the conformal gauge. 
The above two equations reduce to the corresponding 
continuity and Euler equations as shown in our previous 
work ^ when the entropy perturbation vanishes. For 
the dark fluid in this paper, we assume the shear pertur- 
bation ffd — 0. In our calculations, the adiabatic initial 
conditions will be taken. To perform the numerical cal- 
culation, we modified the publicly available cosmoMC 
package 12411 to include the dark fluid perturbation in the 
CAMB code which is used to calculate the theoret- 
ical CMB power spectrum. We added new parameters 
and modified the perturbations equations for the dark 
energy in the CAMB code. 



III. CONSTRAINT METHOD AND RESULTS 

A. Implications on CMB anisotropy for model 
parameter (?s,eff 

In our previous paper Q , for the only adiabatic pertur- 
bation case, we have discussed the implication on CMB 
anisotropic power spectra of the model parameter a and 
Bs, for the details please see [1]. When we fix the val- 
ues of the adiabatic and efi^ective sound speeds, the other 
relevant model parameters have the same effects to the 
CMB power spectra as shown in our previous paper Q. 
So, in this paper, we focus on the model parameter e//- 
To illustrate how the CMB temperature anisotropics are 
characterized by different values of gj^, we choose dif- 
ferent values of in the range [0, 1] with the other 
relevant cosmological models fixed to the mean values 
obtained in 



For different values of the effective sound speed in 
the range [0,1], we plotted their effects on the CMB 
anisotropic power spectra in Figure [TJ As shown in this 
figure, CMB power spectra favor small values of cleft- 



(11) 



(12) 

I 




FIG. 1. The effects on CMB temperature power spectra 
of model parameter c^.^ff- The solid red, blue dashed, 
green dotted and yellow dash-dotted lines are for c^^^ff = 
0,0.1, 0.5, 1 respectively, where the other relevant cosmologi- 
cal parameters are fixed to their mean values obtained in Ref. 
[3. The CMB power spectrum shows that small values of 
Cs,e// are favored. 

For large values of ^jj: , the gravitational potential de- 
cays fastly due to pressure support of the UDF fluctua- 
tions during UDF domination. The effect monotonically 
decrease as decrease to zero which was analyzed in 

Ref. We' expect the CMB power spectra can give a 
tight constraint to the model parameter . 

B. Method and data points 

To constrain the model parameter space, we use the 
Markov Chain Monte Carlo (MCMC) method. We 
adopted the following 8-dimensional parameter space 

P = {iUb, Bs, T, a, B„ c2 g^^, n„ log[10i%]}. (13) 

The priors of the model parameters are shown in Table 
m We adopted fc^o = 0.05Mpc~^ as the pivot scale of the 
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initial scalar power spectrum. 



Model Prameters 


Priors 




[0.005, 1] 


es 


[0.5, 10] 


T 


[0.01,0.8] 


a 


[-1,1] 




[0,1] 




[0,1] 


ris 


[0.5, 1.5] 


log[10^°A] 


[2.7,4] 


Age 


lOGyr <to< 20Gyr 




0.022 ± 0.002 [26] 


Ho 


74.2 ± 3.6kms-^Mpc-^ [27] 



TABLE I. The priors for the cosmological model parameters 
and other priors. The pivot scale of the initial scalar power 
spectrum kso = 0.05Mpc~^ is used in this paper. 



The total likelihood jC (x e ^ 1"^ should be calculated 
to get the model parameter space, where is given as 



2 _ 2 I 2 

X — XCMB + Xbao 



XSN- 



(14) 



For CMB data set, the temperature power spectrum from 
WMAP 7-year data [H] are emplwed. For the BAO 
information, the SDSS data points [23| are used. For SN 
la, we use the 580 Union2.1 data sets with systematic 
errors jl3]- For the detailed description, please see Refs. 

aim. 



C. Fitting Results and discussion 

We ran 8 chains in parallel on the Computational Clus- 
ter for Cosmos (3C) and checked the convergence {R—l 
is of the order 0.01). The obtained results are shown in 
Table inland Figure H 




FIG. 2. The 1 — D likelihood and 2 — D contour plots for the model parameters with 1, 2, 3(t regions. 



r 



The most relevant and interesting model parameters c^s ad ~ ^ effective sound speed c^e//- They 

to the UDF are the constant adiabatic sound speed 
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Prameters 



Mean with errors 



Qs 

T 

a 

Us 
fid 

Age/GyT 
fib 

Zre 

Ho 



0.0228 



+0.000628+0.00128+0.0020315 
-0.000628-0.00120-0.00176 



1 n/IQl^+0-00274+0. 00553+0. 00844 
l.U^»O_QQQ27i_0. 00522-0. 00780 
n nQ/lQ+0 00690+0. 0282+0. 0522 
U.Ut>y:a_o. 00837-0. 0253-0. 0418 
n nnnQc;7+0-0018''+0-00345+0. 00523 
U.UUUyO I -0,00184-0.00406-0.00675 
n 917+0.0167+0.0359+0.0573 
'J-^J^ ' -0.0168-0.0308-0.0457 

n nni cc+o. 000319+0.00241+0. 00493 

U.UU100_o 00155-0.00155-0.00155 
n QQQ+0. 0230+0. 0509+0. 0811 
U.yyy_o. 0227-0. 0393-0. 0548 
o nyni +0.0370+0.0730+0.111 
O.U I ^1-0.0362-0.0713-0.112 
n Qce+0. 00249+0. 004673+0. 00679 
U.yaO_o. 00252-0. 00530-0. 00840 
10 790+0.147+0.295+0.446 
10. / ^0-0.145-0.280-0.435 
n n/l'RQ+0-00252+0. 00530+0. 00841 
U.U10»_Q 00249-0.00467-0.00679 
in QO/l +1.260+2.498+3.873 
lU.yzy:_i 242-2.410-3.824 
70 1 OQ+1. 824+3. 626+5. 537 
' ^-^"^-1.816-3.631-5.311 



TABLE II. The cosmological model and derived model 
parameters space with 1,2,3a regions where SN 
Union2.1+BAO+CMB data sets are used. 




FIG. 3. The evolutions of the EoS Wd{a) for UDF with respect 
to the scale factor a, where the mean values of the relevant 
model parameters are adopted. 



equal to each other when the entropy perturbation van- 
ishes. The results show that the currently available data 
sets from SN, CMB and BAO can give a tight constraint 
to the model parameter space and favor a model with 

r2 - n nm ^^+0.000319+0.00241+0.00493 1 o 

'^s,eff — ^•^'-'^"J'J-0. 00155-0. 00155-0. 00155 1" i, oa le- 

gions respectively. When the entropy perturbation is in- 
cluded, a negative adiabatic sound speed is favored which 
is different from that of the pure adiabatic case. By us- 
ing the obtained mean values, we plotted the evolutions 
of EoS for the UDF with respect to the scale factor a 
in Figure |31 From the EoS figure, one can read off that 
the UDF behaves like cold dark matter at the earlier 
epoch and like dark energy at the later epoch. Then the 
small perturbations of UDF can grow into the large scale 
structure of our Universe. We also show the Cf^ power 
spectra for ACDM model and observed data points in 
Figure |4] where the mean values of relevant model pa- 
rameters are adopted. It implies that current cosmic 
observational data points can not discriminate ACDM 
model from the UDF model. Though the data points 
favor a small amount of entropy perturbations, but the 
pure adiabatic perturbations are not rule out. 




















H 
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FIG. 4. The CMB Cf^ power spectrum v.s. multiple moment 
/, where the grey dots with error bars denote the observed 
data with their corresponding uncertainties from WMAP 7- 
year results, the red dashed lines are for the UDF model 
with mean values as shown in Table [ill the blue solid lines 
are for ACDM model with mean values taken from [28] with 
WMAP-fBAO-}-i?o constraint results. The bottom panel 
shows the ratios to ACDM model. 



We also show the evolutions of 5d and Vd of the 
UDF with respect to the redshift z on the scale k — 
10~'^Mpc~^ in Figure [5l One can easily see the growth 
of the UDF perturbations with the evolution of our Uni- 
verse. 



IV. SUMMARY 



In this paper, we continued to study a unified dark 
fluid model with a constant adiabatic sound speed but 
with the entropic perturbations. As comparisons to the 
previous work we took the effective sound speed as 
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FIG. 5. Evolution of 5d (thin) and Vd (dashed) with respect 
to the redshift z where the mean values of model parameters 
listed in Table |ll] are adopted. Here k = 10"^Mpc"^. 

a free model parameter to characterize the micro-scale 
properties in additional to its EoS instead of fixing it 
to a specific value by hand. And then we try to use the 



currently available cosmic observations which include SN 
Union2.1, BAO and full information of CMB to deter- 
mine the properties of this UDF model via the MCMC 
method. A tight constraint was obtained as shown in 
Table HT] The results show that a UDF model with small 
entropy perturbation is favored but the pure adiabatic 
case is not ruled out. We also show the evolutions of the 
EoS of the UDF and its perturbations with the evolu- 
tion of our Universe. The analysis tell us that the UDF 
behaves like cold dark matter at early epoch and like 
dark energy at late time. One can also see the growth of 
the perturbations of the UDF with the expansion of our 
Universe in Figure [5] It shows the possibility for a small 
UDF perturbation to grow into a large scale structure of 
the Universe. We expect our study can shed light on the 
understanding of the dark side of our Universe. 
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